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SECTION  I 


INTRODUCTION 

The  expanding  utilization  of  laminated  composite  materials  in  air- 
craft structural  components  has  led  engineers  to  seek  new  ways  to  utilize 
fully  these  light-weight  materials.  One  aspect  of  this  search  has  been 
the  study  of  aeroelastic  tailoring  of  composite  materials.  Aeroelastic 
tailoring  of  a lifting  surface  involves  the  synthesis  of  a composite  ma- 
terial structure  whose  special  directional  properties  enhance  the  aero- 
elastic performance  of  the  airfoil.  Primary  attention  has  been  focused 
upon  improving  conventional  aircraft  designs  such  as  sweptback  wing  con- 
figurations. The  purpose  of  this  report  is  to  illustrate  the  potential 
benefits  to  be  gained  by  applying  the  aeroelastic  tailoring  concept  to  a 
nonconventional  design,  an  aircraft  with  a forward  swept  wing. 

A report  by  Diederich  and  Budiansky  [1]  discusses  the  aeroelastic 
divergence  of  metallic  swept  wings.  Their  report  shows  convincingly  that, 
from  a static  aeroelastic  stability  point  of  view,  sweptback  metallic  wings 
are  far  superior  to  sweptforward  metallic  wings.  As  a result  of  their 
study,  the  forward  swept  metallic  wing  has  not  been  a serious  design  con- 
cept for  some  time. 

Although  the  metallic  forward  swept  wing  is  inferior,  from  a structural 
stiffness  standpoint,  the  forward  swept  wing  possesses  some  aerodynamic 
advantages  over  the  sweptback  wing.  Knight  and  Noyes  [2]  noted  in  1931 
that  wings  with  20°  of  forward  sweep  experience  stall  at  higher  angles  of 
attack  than  similar  20°  sweptback  wings.  Thus,  sweptforward  wings  may 
develop  higher  lift  coefficients  tnan  similar  sweptback  wings.  In  addition, 
Jones  [3]  notes  that  sweepback  may  produce  a pronounced  rolling  or  pitching 
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instability  because  of  the  premature  tip  stall  of  these  wings.  Con- 
versely, at  high  angles  of  attack,  forward  swept  wings  begin  to  stall  first 
at  the  wing  root  and  are  not  subject  to  this  phenomenon.  On  the  other 
hand,  the  inherent  lack  of  directional  stability  of  a forward  swept  fly- 
ing wing  makes  it  necessary  to  provide  a generous  vertical  stabilizer  and 
rudder.  Little  has  been  published  concerning  the  stability  and  control 
characteristics  of  forward  swept  wing-fuselage  configurations. 

The  first  aircraft  to  employ  a forward  swept  wing  was  the  Junkers 
Ju  287  jet  bomber  [4]  developed  in  Germany  during  World  War  II.  This  test 
aircraft  flew  on  16  occasions  in  1944,  developing  speeds  of  over  650  km/hr 
(404  mph).  The  wing  on  the  Ju  287  was  swept  forward  nearly  25°,  primarily 
to  provide  good  low  speed  flight  characteristics.  Limited  flight  tests, 
brought  to  a halt  by  advancing  Soviet  armies,  did  not  disclose  any  severe 
aeroelastic  or  directional  stability  problems.  The  interested  reader  will 
find  the  information  in  Reference  4 a fascinating  account  of  a little 
known  chapter  in  aircraft  history. 

A modern-day  example  of  an  airplane  with  a forward  swept  wing  is  pro- 
vided by  the  German  Hansa  twin- jet  executive  transport  [5].  This  aircraft 
is  the  only  forward  swept  wing  aircraft  other  than  the  Ju  287  that  has 
actually  flown.  Developed  by  Messerschmitt-Bolkow-Blohm,  the  prototype 
HFB  320  Hansa  first  flew  in  1964.  The  metallic  wing  on  the  HFB  320  ex- 
hibits an  aspect  ratio  of  6,  a taper  ratio  of  1:3,  and  employs  15°  of  for- 
ward sweep  at  the  quarter-chord  line  to  avoid  the  necessity  of  taking  the 
main  spar  of  the  wing  through  the  passenger  cabin.  The  aircraft  has  a max- 
imum cruise  speed  of  825  km/hr  (513  mph)  at  25,000  ft.  A stretched  ver- 
sion of  the  HFB  320,  designated  as  the  HFB  330  Hansa  has  also  been  developed. 
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Renewed  ircerest  in  forward  swept  wing  technology  is  due  to  N.  J. 

Krone,  Jr.  [0].  Reference  6 illustrates  the  feasibility  of  the  forward 
swept  wing  concept  when  laminated  composite  materials,  rather  than  metals, 
are  used  in  the  wing  structure.  Krone's  conclusion,  arrived  at  through 
a comprehensive  set  of  numerical  examples,  is  that  the  use  of  laminated 
composite  materials  in  the  wing  structure  offers  the  design  flexibility 
to  overcome  the  weight  disadvantages  encountered  in  conventional  metallic, 
forward  swept  wing  design.  Thus,  proper  aeroelastic  tailoring  appears 
to  lessen  the  severity  of  the  aeroelastic  divergence  problem  of  forward 
swept  wings. 

The  results  in  Reference  6 leave  unsettled  several  questions  about 
the  basic  mechanism  by  which  static  aeroelastic  stability  improvement  is 
achieved  through  the  use  of  laminated  composites.  The  present  report  pro- 
vides a study  of  the  static  aeroelastic  divergence  characteristics  of  wings 
whose  structure  is  constructed  solely  of  laminated,  fibrous  composite  materi- 
als. This  study  uses  a mathematical  idealization  of  the  wing  divergence 
problem  similar  to  that  used  in  Reference  1,  but  incorporates  laminated 
beam  theory  into  the  analysis.  Attention  is  focused  upon  the  additional 
coupling  between  bending  and  torsional  deformations  of  a wing,  introduced 
by  the  laminated  composite  construction,  and  its  influence  on  swept  wing 
divergence. 

The  mathematical  model  used  throughout  the  study  employs  a strip 
theory  representation  of  the  aerodynamic  loads  that  arise  from  wing  bending 
and  torsional  deformation  of  high  aspect  ratio  wings.  This  representation, 
together  with  laminated  beam  theory,  allows  the  solution  of  the  wing  diver- 
gence problem  in  analytical  form.  Such  a solution  is  useful,  since  the 
effect  of  geometrical,  structural  and  aeroelastic  parameters  can  be  seen 
readily. 
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SECTION  II 

AEROELASTIC  MODEL  IDEALIZATION 

Son.3  insight  into  the  influence  of  laminated  composite  materials 
upon  swept  wing  design  can  be  found  through  a mathematical  analysis  of 
an  idealized  model  of  a slender,  high-aspect-ratio  wing,  such  as  that 
shown  in  Figure  1.  This  section  describes  the  assumptions  used  to  model 
both  the  structural  deformation  behavior  of  the  wing  and  the  aerodynamic 
loads  due  to  this  deformation.  This  idealized  model  will  be  used  to  devel- 
op closed-form  solutions  for  aeroelastic  divergence  of  swept,  composite 
wings. 


Figure  1 - Geometry  and  sign  convention  used 
for  the  constant  chord  swept  wing 
analysis. 
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1 . WING  STRUCTURAL  DEFORMATION  MODEL 

The  wing  shown  in  Figure  i is  assumed  to  derive  all  of  its  bending 
and  torsional  stiffness  from  thin,  laminated  composite,  cover-sheets  that 
form  the  upper  and  lower  surfaces  of  the  wing.  It  is  assumed  that  the 
resulting  box-beam  wing  structure  is  such  that  its  deformation  may  be  repre- 
sented by  a bending  deflection  h(y),  positive  upward,  along  a straight 
reference  axis  (the  y-axis  in  Figure  1)  and  a rotation  a(y),  positive  nose- 
up,  about  this  axis.  In  addition,  it  is  further  assumed  that  wing  chord- 
wise  sections,  perpendicular  to  the  beam's  reference  axis,  are  rigid  so 
that  wing  deformation  is  a function  only  of  the  spanwise  coordinate,  y. 

The  enforcement  of  this  type  deformation  is  essential  to  the  work  that 
follows.  A similar  structural  model  is  used  by  Housner  and  Stein  [7]  to 
study  the  flutter  of  sweptback  laminated  composite  wings. 

Using  the  box-beam  model,  equivalent  bending  and  torsional  stiffnesses 
of  the  wing,  El  and  GJ,  respectively,  may  be  computed  employing  the  clas- 
sical Euler-Bernoulli  assumptions.  As  is  suggested  in  Reference  7,  the 
derivation  of  the  equivalent  El  and  GO  for  the  structure  is  best  accomplished 
through  strain  energy  methods.  The  Euler-Bernoulli  displacement  behavior 
constrains  the  individual  laminae  in  the  upper  and  lower  laminated  cover 
sheets  to  act  as  a unit  and  supplies  the  strain-displacement  relations  for 
the  wing  structure.  The  definition  of  the  stress-strain  relations  of  each 
of  the  laminae  furnishes  the  additional  information  necessary  to  formulate 
the  strain  energy  expression  for  the  box-beam. 

Attention  will  now  turn  to  a summary  of  the  development  of  the  ex- 
pressions for  El  and  GJ.  Since  terms  that  involve  the  individual  lamina 
stress-strain  properties  appear  in  the  expressions  for  El  and  GJ,  the 
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development  of  the  stress-strain  relations  for  the  lamina  will  be  reviewed. 
A full  discussion  of  the  development  of  the  equilibrium  equations  that 
govern  the  laminated  beam  is  contained  in  Appendices  A and  B. 


Figure  2 - Composite  material  lamina  principal 
axes,  1 and  2,  with  respect  to  the 
wing  reference  axes,  x and  y. 

Each  lamina  of  the  box-beam  is  itself  orthotropic  with  respect  to  a 
set  of  principal  axes,  one  axis  being  oriented  in  the  direction  of  the 
lamina  fibers,  while  the  other  axis  is  perpendicular  to  these  fibers.  These 
principal  axes,  denoted  as  axis  1 and  axis  2 and  shown  in  Figure  2,  lie  in 
a plane  parallel  to  the  x-y  reference  plane  of  the  wing,  but  are  oriented 
at  an  angle  e with  respect  to  the  x-y  system.  The  angle  e is  measured 
positive  in  the  counterclockwise  direction  from  the  x-axis  shown  in  Figure  2. 


As  shown  in  Reference  8,  the  relationship  between  inplane  stresses 
and  strains  may  be  written,  for  a single  orthotropic  lamina,  as 


6 


The  terms  Q^.  are  functions  of  e,  the  fiber  angle,  and  the  orthotropic 
engineering  constants  Q^.  given  by  the  following  expressions. 


Qn  = ' wi2  ^21  ^ 

(2a) 

G12  = u12  " y12  u21^ 

(2b)’ 

= y2i  ~ ^21  u12^ 

Q22  s E2/(l  - p12  y21) 

(2c) 

Q66  = G12 

(2d) 

The  terms  Q. . are  defined  in  a convenient  manner  by  first  defining  a 

* J 

set  of  parameters  U-j  that  are  functions  of  Q^-  [8]. 

U1  = [3Qn  + 3Q22  + 2Q1Z  + 4Q66]/8 

(3a) 

U2  = W11  " W* 

(3b) 

U3  * Kn  + 022  ' 4IW/8 

(3c) 

U4  ’ [Q11  + ^22  + 8^12  " 

(3d) 

U5  ' W-,1  + Q22  ' ^12  + 4W/8 

(3e) 

In  terms  of  these  new  parameters  U.  and  the  angle  e,  the  (L. 

1 • vJ 

1 terms 

become 

^11  = U1  + U2cos2e  + U3cos40 

(4a) 

^12  = U4  " U3cos4e 

(4b) 

^22  = U1  " U2cos2e  + U3cos4e 

(4c) 

^16  = I U2sin2e  + U3S™40 

(4d) 

^26  = \ U2s^n2e  " U3sin4e 

(4e) 

^66  = U5  ‘ U3cos4e 

(4f) 
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Since  the  individual  laminae  are  constrained  to  act  as  a unit,  in  a 
manner  prescribed  by  the  Euler-Bernoulli  deformation  assumptions,  a strain 
energy  functional  can  be  constructed  by  summing  the  individual  strain 
energies  of  each  lamina.  The  reference  surface  for  this  strain  energy  ex- 
pression is  the  middle  surface  of  the  wing  box.  The  total  strain  energy, 

denoted  as  U,  is  written  as  follows  (note  that  ( )'  = d{  )/dy): 

£ 

“■  / ?EIo<h">2  - Vh'v)  + ?W>2 

0 j (5) 

*b33<vov>  - B2z<h"vo'> + i sj0(“'>z{dy 

The  terms  multiplying  the  deflection  derivatives  in  Equation  (5)  are  defined 
in  Appendix  A.  These  terms  are  linear  functions  of  the  wing  chord  and  also 
involve  the  Q^.  parameters  of  the  individual  laminae. 

From  Equation  (5)  it  is  seen  that  the  bending  deformation  h(y)  is  elas- 
tically coupled  both  to  the  torsional  deformation  a(y)  and  to  the  exten- 
sional  deformation  of  the  middle  surface,  v0(y).  The  extensional  deforma- 
tion v (y)  will  be  zero  if  the  middle  surface  is  also  the  neutral  surface 
o 

for  bending,  a situation  that  occurs  if  the  upper  and  lower  laminated  cover 
sheets  are  oriented  symmetrically  with  respect  to  the  middle  surface. 

The  Principle  of  Virtual  Work  is  used  to  derive  the  differential 
equations  governing  the  static  equilibrium  of  the  laminated  beam  under 
the  action  of  the  upward  distributed  load  p(y)  and  a distributed  torque 
t(y).  These  equations,  derived  in  Appendix  B,  are 

[Elh"  - KoT  = p(y)  (6) 

[-KhM  + GJa']1  = -t(y)  (7) 

Note  that  in  Equations  (6)  and  (7)  ( )'  * d(  )/dy.  The  terms  El,  GJ  and  K 


-£^r, v7: 73%  ^ ^ 


are  related  to  the  parameters  in  Equation  (5)  by  the  expressions 


El  = EI0-(B22)VA22 

(8a) 

GJ  = GJ0-(B33)2/A22 

(8b) 

K = Ko“B22B33/,A22 

(8c) 

Equations  (8a,b,c)  are  the  result  of  the  elimination  of  the  vQ  degree  of 
freedom  from  the  problem,  an  operation  described  in  detail  in  Appendix  B. 

Equations  (6)  and  (7)  are  coupled  together  by  the  parameter  K,  a 
term  that,  in  general,  is  nonzero  for  a laminated  material.  The  K parameter 
provides  an  elastic  coupling  between  bending  and  torsion  not  found  in  metal- 
lic beam  analysis.  Since  K depends  upon  both  the  stacking  sequence  used 
in  constructing  the  laminate  and  the  lamina  fiber  orientation  angles,  an 
opportunity  to  tailor  the  wing  structure  to  provide  advantageous  bending- 
torsion  coupling  is  presented. 

Boundary  conditions  for  the  cantilevered  wing  consist  of  the  displace- 
ment and  slope  conditions  at  the  wing  root, 

o(0)  = h(0)  = ^ (0)  = 0 (9a,b,c) 

and  the  bending  moment,  shear  and  torque  conditions  at  the  tip,  written  as: 


M(*)  = Elh"  - Ka'  = 0 

(10a) 

VU)  = [Elh”  - Ka']'  = 0 

(10b) 

TU)  = GOa'  - Kh"  = 0 

(10c) 

2.  IDEALIZATION  OF  WING  AERODYNAMIC  LOADS 

Aerodynamic  strip  theory  provides  a meaningful  and  somewhat  accurate 
way  of  describing,  in  mathematical  terms,  the  air  loads  acting  upon  a 
slender  wing  as  the  result  of  structural  deformation.  A complete  discussion 
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of  this  aerodynamic  formulation  is  given  in  Reference  9. 

The  aerodynamic  forces  and  moments  acting  upon  wing  segments  that  are 
perpendicular  to  the  reference  axis  (the  y-axis)  of  the  wing  may  be  com- 
puted, provided  an  effective  wing  root  and  tip  are  assumed,  an  assumption 
that  may  involve  some  error.  However,  the  strip  theory  approach  produces 
simple  equations  that  have  been  shown  to  give  satisfactory  results  for 
large  aspect  ratio  metallic  wings. 


EFFECTIVE  "V 
TIP 


Figure  3 - Slender  swept  wing  geometry  used  to 
describe  aerodynamic  loads,  showing 
chordwise  segments  perpendicular  to 
the  swept  reference  axis. 
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To  formulate  the  divergence  problem,  the  only  aerodynamic  loads  to  be 
considered  are  those  caused  by  deformation  h(y)  and  a(y).  Figure  3 defines 
tt’e  geometrical  parameters  used  to  define  these  aerodynamic  loads.  In 
terms  of  the  upward  deflection  of  the  reference  axis,  h(y),  and  the  tor- 
sional rotation,  a(y),  of  the  wing  chordwise  sections,  about  this  axis,  the 
lift  force  per  unit  length,  positive  upward,  is 


O JL 

P(y)  = qcaQcos  a(q  - tan  a) 


(ID 


The  torque  per  unit  length  due  to  structural  deformation  is  expressed  as 


t(y)  = qceaQcos2A(a  - ^ tan  a)  (12) 

The  term  aQ  appearing  in  Equations  (11)  and  (12)  is  the  two-dimensional 
lift  curve  slope  of  the  wing  sections. 

The  lift  curve  slope  can  be  empirically  corrected  to  reflect  the 
effects  of  finite  span,  wing  sweep  and  compressibility.  The  present  report 
makes  no  modification  of  the  2-D  lift  curve  slope  to  account  for  compres- 
sibility. However,  the  2-D  lift  curve  slope  is  empirically  corrected  to 
take  into  account  aspect  ratio,  AR,  and  sweep  angle  A by  applying  the  cor- 
rection 


ao  “ ao(A=0) 


AR 


(13) 


AR  + 4cosA 

Equation  (13)  is  adapted  from  that  used  in  Reference  1. 

Equations  (11)  and  (12)  are  inserted  into  Equations  (6)  and  (7),  re- 
spectively, to  provide  the  differential  equations  governing  the  aeroelastic 
divergence  of  composite  swept  wings.  The  solution  to  these  equations  for 
several  cases  is  discussed  in  the  ensuing  sections. 


3.  DIVERGENCE  OF  A UNIFORM  CHORD  WING  WITH  BENDING  FLEXIBILITY  ONLY 

The  simplest  wing  divergence  case  to  consider  is  one  in  which  the  wing 


S^V?-  ^j5! e35^ 


has  constant  geometrical  and  stiffness  properties  along  its  span  and  has 
no  torsional  flexibility.  For  this  case,  the  governing  equation  of  equilib- 


rium becomes 


£Ihiv  + qcaQ(^)  sinACOSA  = 0 


Equation  (14)  can  be  transformed  into  the  following  nondimensional  form: 


- br  = 0 


by  defining  the  new  variables 


5 = 1 - y/, 


b = qcao£3sinAcosA 


Boundary  conditions  for  this  cantilever  wing  problem,  in  terms  of  the  vari- 


able 5,  are: 


r(D  = 0,  MO I - o <rr(p)  _ 0 
d«  ’ dE2 


(16a, b,c) 


This  problem  is  an  eigenvalue  problem  and  is  discussed  at  length  in 
Reference  10  (pp.  311  - 317).  The  solution  to  this  problem  yields  a critical 
value  of  dynamic  pressure,  q,  at  divergence,  denoted  as  qQ,  given  by  the 


following  expression: 


6.33  El 

3 

a a |sinA|cosA 


Equation  {17}  is  valid  only  if  A,  the  sweep  angle,  is  less  than  zero;  this  value 
of  a occurs  when  the  wing  is  sweptforward.  A sweptback  wing  (a>0)  will  not 
diver.-  .nder  toe  present  assumptions. 

With  no  torsional  deformation  permitted,  the  dynamic  pressure  at  diver- 
gence is  maximized  by  maximizing  the  parameter  El.  The  maximum  value  of 


saa; 


the  parameter  El  is  attained,  in  the  case  of  a fibrous  composite,  by 
aligning  the  fibers  parallel  to  the  swept  reference  axis.  This  result  is 
not  surprising  and,  indeed,  is  somewhat  discouraging,  because  the  solution 
to  designing  such  a structure  lies  In  adding  more  composite  material  as 
the  wing  is  sweptforw^rd . This  strategy  is  identical  to  that  followed  for 
conventional  metallic  designs. 

4.  DIVERGENCE  OF  A SWEPT  WING  WITH  UNIFORM  GEOMETRICAL  AND  STRUCTURAL 
PROPERTIES  - BENDING  AND  TORSION  INCLUDED 
a.  Equilibrium  Equations 

The  coupled  differential  equations  that  govern  aeroelastic  divergence 
of  a constant  chord,  laminated  composite,  swept  wing  with  both  bending  and 
torsional  flexibility  are  written  as: 

EIh!V  - Ka"1  -1  (qca0cos2A)  (a  - h'tanA)  (18) 

and 

GJa"  - Kh"'  = -(qca  e cos2a)  (a-  h'tan/i)  (19) 

where  ( )'  ~ d(  )/dy. 

To  nondimensionalize  the  above  equations  and  thus  facilitate  the  solu- 
tion, three  new  variables  are  defined: 

n = y/0,  w = h/  , r = dh  = dw  (20a,b,c) 

1 1 dy  dn 

Tnese  new  variables  permit  Equations  (18)  and  (19)  to  be  expressed  as: 

r'"  - k a'*'  = a(a  - r tanA)  (21) 

a"  - gr"  - -b(a  - r tanA)  (22) 

where  now  ( )'  - d(  )/dq.  Parameters  in  the  above  expressions  are  defined 
as  follows: 
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a = qcs,  aQcos  a/EI 


2 2 

b = gee*  aQcos  a/GJ 


(23a) 
(23b) 

U 

k = K/EI  (23c) 

g = K/GJ  (23d) 

The  boundary  conditions  for  these  equations  are  written  as  follows: 


r(0)  = 0 

a(0)  = 0 

r * ( 1 ) - ka'(l)  = 0 
a'(l)  - gr*(l)  = 0 

l"(l)  - ka"(l ) = 0 


(24a) 

(24b) 

(24c) 

(24d) 

(24e) 


Equations  (21)  and  (22)  can  be  combined  into  a single  equation  by  de- 

(25) 


fining  a new  variable,  afi,  given  by  the  relationship 


a = a - T tanA 
e 

The  variable  aQ  represents  the  angle  of  attack  of  a chordwise  section  in- 
duced by  the  torsional  and  bending  deformations  of  the  wing.  The  steps 
necessary  to  accomplish  the  combination  of  these  two  equations  into  a single 
equation  are  outlined  in  Appendix  C.  The  resulting  equation  is  a third 
order  differential  equation  with  constant  coefficients. 


a. 


b(l-k  tanA) 

~ 1 - 

a(g  - tanA) 

0 - kg) 

e 

(l  “ kg) 

°e  = 0 


Equation  (26)  has  boundary  conditions  given  by 

«e(0)  = 0,  a;(l)  = 0 


and 


(26) 

(27a, b) 
(27c) 


°e"(1)  + Ck(V'A9~1]  “e  (1)  * 0 

If  the  parameters  k and  g are  zero,  Equation  (26)  and  its  associated  bound- 
ary conditions  reduce  to  those  developed  in  Reference  1 for  the  study  of 
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divergence  of  a metallic  wing. 

Two  aeroelastic  parameters  appear  in  Equation  (26).  These  parameters 
are  defined,  for  ease  of  notation,  as 


a = 


d = 


1-k  tanA 

1 - kg 

tanA  - q 

1 - kg 

_ 

qce£2aQcos2A  1 


(28) 

(29) 


b.  Solution  to  the  Divergence  Equations 

The  solution  to  the  eigenvalue  problem  defined  by  Equation  (26)  and 
its  boundary  conditions.  Equations  (27a, b,c), follows  that  given  in  Refer- 
ence 1.  This  solution  is  summarized  in  Appendix  C of  the  present  report. 

Figure  4 displays  the  interdependence  of  the  critical  values  of  a 
and  d found  for  this  problem.  The  region  above  the  solid  curve  shown  in 
Figure  4 represents  a region  in  which  combinations  of  values  of  a and  d 
result  in  a divergence  condition  for  the  swept  wing.  Values  of  a and  d 
below  the  solid  curve  correspond  to  those  values  of  the  aeroelastic  param- 
eters that  are  associated  with  a stable  wing.  Thus,  the  solid  curve  cor- 
responds to  a stability  boundary  for  the  constant  chord  wing. 

As  discussed  in  Reference  1,  an  approximation  to  the  curve  shown  in 
Figure  4 provides  an  accurate  approximate  relation  between  a^  and  d^,  the 
values  of  a and  d at  divergence.  This  approximation  is  developed  by  first 
recognizing  that  a is  equal  to  zero  if  the  wing  is  rigid  in  torsion.  This 
condition  corresponds  to  bending  divergence  of  a swept  wing,  discussed  in 
Section  II,  part  3.  In  this  case, 


dQ  = -J-  = -6.333 


(30) 


If  the  wing  is  rigid  in  bending,  but  has  torsional  flexibility,  then  d = 0 


15 


and  the  exact  solution  to  the  divergence  problem  results  in  a value  of  a^ 
given  by 

aD  = (,/2)2  (31) 

If  a linear  algebraic  relationship  between  ?p  and  dp  is  assumed,  taking 
into  account  the  results  in  Equations  (30)  and  (31),  an  approximate  relation- 
ship between  ap  and  dp  can  be  developed.  This  relationship  is  given  by  the 
equation 


The  approximate  relationship  given  by  Equation  (32)  is  shown  as  a 
dashed  line  in  Figure  4.  For  the  range  of  values  shown  in  this  figure,  the 
agreement  between  the  exact  solution  and  the  approximate  solution  is  excel- 
lent. 


Combining  Equations  (28)  and  (29)  with  Equation  (32),  the  dynamic 
pressure  at  divergence  is  found  to  be 


• 

2.47  (1  - kq) 

El 

o 2 

ciTa0cos  a 

1 - k tanA  - 0.39(tanA  - g) 
_ (?)  # 

Equation  (33)  provides  useful  information  about  the  divergence  be- 
havior of  the  swept  wing.  First  of  all,  if  qp  is  a negative  number,  no  real 
divergence  speed  occurs,  since  the  airspeed  V must  be  an  imaginary  number 
if  qp  is  to  be  less  than  zero.  Therefore,  combinations  of  parameters  that 
result  in  negative  values  of  qp  are  desirable.  From  strain  energy  consid- 
erations, the  factor  (1  - kg)  in  Equation  (33)  will  always  be  a positive 
number.  Thus,  the  denominator,  containing  terms  dependent  upon  A and  the 
geometric  and  structural  characteristics  of  the  wing,  controls  the  sign  of 


q 


D* 
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The  divergence  dynamic  pressure  q^  changes  sign  at  infinity.  The 
parameter  q^  becomes  infinite  when  the  following  relationship  is  satisfied: 


tan  a 


cr 


, + 2'56(?) 


(34) 


/ 

Values  of  tan  a larger  than  that  given  in  Equation  (34)  result  in  finite  / 

f 

/ 

negative  values  of  qQ  that,  in  turn,  yield  imaginary  divergence  velocities. 

The  value  of  a determined  in  Equation  (34)  is  denoted  as  a„„  because  it 

cr  / 

represents  that  value  of  the  sweep  angle  above  which  swept  wing  div/rgence 

/ 

is  theoretically  impossible.  / 

To  interpret  and  to  compare  the  information  provided  by  Equation  (34), 

/ 

it  is  noted  that  a conventional  metallic  wing  has  the  parameter  g equal  to 

zero.  Thus,  a i's  a constant  for  the  metallic  wing.  For/a  laminated  com- 
er 

posite  wing,  Acr  will  be  a function  of  the  laminate  characteristics.  Fur- 
thermore, for  a metallic  wing  operating  at  subsonic  SDeeds,  a will  be 

V I 

greater  than  zero  for  the  usual  case  in  which  the  parameter  e is  greater 
than  zero.  This  means  that  the  metallic  wing  must  be  swept  back  to  pre- 
clude divergence.  On  the  other  hand,  with  a composite  material  wing  for 
which  the  parameter  g is  less  than  zero,  the  possibility  exists  that  A 
is  less  than  zero.  In  chis  case,  sweep  angles  exist,  in  the  forward  swept 
range,  for  which  divergence  will  not  occur  at  any  flight  speed, 
c.  Example  Cases 

An  example  of  the  use  of  Equation  (34)  to  study  divergence  is  provided 


by  two  related  wing  configurations.  In  the  first  case,  all  of  the  laminate 
fibers  of  the  cover  sheets  are  oriented  along  a common  direction,  denoted 
as  the  angle  e.  In  the  second  case,  10%  of  the  fibers  are  fixed  at  an 
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angle  e = 0°,  25$  of  the  fibers  are  oriented  symmetrically  at  e = + 45° 
(12.5$  are  + 45°,  12.5$  are  -45°),  while  the  remaining  65$  of  the  laminate 
has  fibers  oriented  at  an  angle  0.  Geometric  and  lamina  material  data  for 
the  wing  are  as  follows: 

Ming  and  Lamina  Properties 

Wing  Aspect  Ratio,  AR  = 6 Offset  between  aerodynamic  center  and 


reference  axis,  e = 0.10c 


Lamina  Material:  Boron-epoxy 
E-|  = 32.5  x 10®psi.  (tension) 
E2  = 3.2  x 1 06psi . (tension) 
G,9  = 1.05  x 106psi. 


= 0.36 


The  ratio  of  box-beair.  depth  to  cover-sheet  thickness  is  20:1. 

Figure  5 shows  the  relationship  between  the  parameter  g=K/GJ  and  0 
for  each  case.  For  Case  1,  in  which  100$  of  the  fibers  can  be  tailored, 
the  absolute  value  of  g has  a maximum  value  of  about  1.20.  Also  to  be 
noted  is  that  g is  an  asymmetrical  function  (about  0 = 90°)  of  0.  In  Case 
2,  where  only  65$  of  the  fibers  are  available  for  design  orientation,  a 
reduction  in  the  percentage  of  fibers  that  can  be  used  for  tailoring  is 
seen  to  result  in  a decrease  in  the  maximum  absolute  value  of  g to  near 
| g | = 0.50.  In  addition,  the  position  of  the  min-max  values  of  g in  the 
latter  case  lies  farther  away  from  the  0 = 90°  position  than  that  found 
for  Case  1.  It  should  be  noted  that  both  El  and  GJ  are  symmetrical  functions 
of  0 (about  0 = 90°). 

Using  Equation  (34),  the  relationship  between  A£r  and  0 for  the  two  pre- 
viously discussed  cases  can  be  determined.  These  results  are  shown  in  Fig- 


ure 6.  In  Case  1,  proper  orientation  of  the  laminate  fibers  can  preclude 
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The  relationship  between  the  elastic  coupling 
parameter  g=K/G0  and  the  fiber  angle  0 for  two 
example  laminates. 


wing  divergence  until  a forward  sweep  angle  of  nearly  49°.  This  occurs  if 
laminate  fibers  are  oriented  at  an  angle  about  12°  ahead  of  the  swept  refer- 
ence axis  (e  = 102°).  For  Case  2,  in  which  fewer  fibers  are  used  in  the 
tailoring  process,  the  results  are  less  dramatic,  but  still  impressive. 

When  e is  about  110°  (20°  forward  of  the  reference  axis)  the  wing  may  still 
be  swept  forward  nearly  28°  before  a real  divergence  speed  is  encountered. 

Figure  7 shows  the  behavior  of  the  divergence  speed  of  the  uniform 
property  wing  in  Case  1 as  the  fiber  angle  changes.  The  results  are 
normalized  with  respect  to  VDg,  the  value  of  VQ  that  results  when  the  wing 
is  unswept  and  all  the  fibers  are  oriented  along  e = 0.  Three  sweep  angles 
are  considered  in  Figure  7;  these  angles  are  A = 0°,  -30°,  -60°,  corres- 
ponding to  an  unswept  wing  and  two  forward  swept  wings.  In  the  range  0°  £ 
e < 90°,  the  divergence  speeds  of  the  wing  for  each  of  the  sweep  angles 
are  relatively  low.  However,  as  8 is  increased  beyond  90°  and  the  fibers 
are  rotated  forward  of  the  wing  reference  axis,  the  divergence  speed  in- 
creases rapidly.  In  the  cases  of  the  unswept  wing  and  the  wing  with  30°  of 
forward  sweep,  the  divergence  speed  becomes  infinite  for  fiber  angles 
e = 95°.  There  is  a range  of  laminate  fiber  angles  for  which  divergence 
is  impossible  in  these  latter  two  sweep  cases.  The  range  of  angles  in 
which  divergence  is  precluded  is  seen  to  decrease  with  increasing  forward 
sweep.  Referring  back  to  Figure  6,  it  is  to  be  expected  that  this  range 
or  bandwidth  approaches  zero  as  the  angle  of  forward  sweep  approaches  49°. 
Returning  to  Figure  7,  it  is  seen  that,  forA=  -60°,  a maximum  divergence 
speed  is  reached  when  e is  approximately  equal  to  100°.  For  this  fiber 
orientation,  this  particular  forward  swept  wing  has  the  same  divergence 
speed  as  an  unswept  wing  with  all  the  fibers  oriented  along  the  reference 
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Figure  6 - The  relationship  between  the  critical  sweep  angle,  A , and  the 
fiber  orientation  angle  e for  two  example  laminates.  Finite 
divergence  speeds  will  exist  for  sweep  angles  that  lie  in  the 
region  above  the  curves.  Uniform  wing. 
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| axis  (e  = 90°). 

\ An  additional,  important  feature  of  Figure  7 is  that  the  divergence 

I behavior  of  the  unswept  wing  is  not  symmetrical  about  e = 90°,  even  though 

: El  and  GJ  are  symmetrical  functions  about  that  position.  For  a metallic 

l wing  at  zero  sweep,  the  divergence  speed  is  independent  of  the  bending  stiff- 

ness, since  bending  deformation  produces  no  lift.  However,  in  the  case  of 
the  laminated  composite  wing,  bending  deformation,  caused  by  lift,  causes 
torsional  deformation  that,  in  turn,  leads  to  changes  in  lift.  This  elastic 
\ coupling  feature  is  clearly  seen  in  Figure  7. 

\ 5.  DIVERGENCE  OF  LINEARLY  TAPERED  WINGS  WITH  GEOMETRICALLY  SIMILAR 

l 

CROSS-SECTIONS 


[ 

; 
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i 
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£ 
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Next  in  order  of  complexity,  after  the  uniform  property  wing,  is  a 
wing  whose  properties  have  a linear  variation  along  the  wing  span.  A solu- 
tion to  Equations  (6)  and  (7)  is  discussed  in  Reference  1 for  the  metallic 
wing  whose  wing  chord  varies  linearly  along  the  span  while  the  bending  and 
torsional  stiffnesses  vary  as  the  fourth  power  of  the  chord.  Such  a case 
occurs  when  a wing  has  geometrically  similar  cross-sections.  For  a com- 
posite wing,  the  parameter  K also  varies  as  the  fourth  power  of  the  chord 
if  both  geometrically  similar  cross-sections  are  used  and  the  chord  varies 
linearly  with  distance  along  the  span. 

With  the  assumptions  of  linearly  varying  chord  and  geometrically  similar 
cross-sections,  expressions  for  the  problem  parameters  are  given  by  the 
following: 

c = fcroot  = fV  e = fer*  EI  = f4EIr’  W = f\' 

(35) 


K = f4Kr 
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where  f = 1 - n(l  - x) 


X = wing  taper  ratio  = ctip/croot 


( ) = parameter  value  evaluated  at  the  wing  root. 
Transforming  Equation  (6)  to  the  independent  variable  f and  then 
nondimensionalizing,  gives  the  bending  equation  (note  that  now  ( )'  = 
d(  )/df): 

f3[rm  - kra,M]  + 8f2[r"  - kr«"]  + 12f[r'  - kr«']  -dTr  + - 0 


where 


v 


qcr*  aosinACOSA 
EIr(l  - X)3 


(36) 


(37) 


Similarly,  the  torsion  equation  is  found  from  Equation  (7)  to  be 
f2[a"  - grr"]  + 4f[a'  - gfr ' ] + aT[a  - rtanA]  = 0 


where 


aT 


qcrer£  aocos  A 
GJr(l  - x)2 


(38) 


(39) 


Boundary  conditions  on  displacement  are,  at  f = 1 (n  = 0), 

a(l ) = 0,  r(l)  = 0 (40a,b) 

The  boundary  conditions  at  f = x (n  = 1)  for  bending  moment,  shear  force 
and  torque  are 


M(x)  = r'(x)  - kra’(x)  = 0 (40c) 

v(x)  = 4x3(r'(x)  - k^'fx))  + x4 ( r" ( x ) -kra"(x))  ■ 0 (40d) 

T(x)  * «'(x)  - gpr'(x)  - 0 (40e) 

Appendix  C,  Section  2,  details  how  Equations  (36)  and  (38)  may  be 
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combined  into  a single  equation  in  terms  of  the  variable  a = a - r tanA. 

e 


This  single  governing  equation  reads  as  follows: 

f 3(5  g + 8f^a“  + (12  + iT)f«;  + (2?t  - dT)ae  = 0 

where  the  aeroelastic  parameters  a^  and  dj.  are  given  by 


aT 


dT 


1 - kftanA 

1 - k gw 

r3r 

L 

tanA  - 

r 

1 ~ kr9r 

qercr^2  a0c°s2A 


2 


GJr(l  - X) 


3 2 

qc^t  a cos  A 
r o 

HIr(l  - X)3 


(41) 


(42) 


(43) 


The  boundary  condition  for  Equation  (41)  at  f = 1 (n  = 0)  is: 

«e(l)  = 0 (44a) 

At  f - x,  the  conditions  that  M(x)  = 0 and  T(X)  = 0 give 

<*;(X)  = 0 (44b) 

The  requirement  that  there  be  zero  shear  at  the  wing  tip,  when  combined 
with  Equation  (38),  evaluated  at  f = X,  provides  the  following  equation. 

X2a"e(x)  + aTae(x)  = 0 (44c) 

Equation  (41)  is  an  Euler  differential  equation  that  can  be  reduced 
to  a linear  differential  equation  with  constant  coefficients  by  changing 
the  variable  f to  f = e*.  The  details  of  solving  Equation  (41)  subject 
to  the  given  boundary  conditions  is  presented  both  in  Reference  1 and  Refer- 
ence 9 (pp.  486  - 487).  These  details  are  omitted  here. 

As  in  the  case  of  the  uniform  property  wing,  the  solution  to  the 
equations  governing  divergence  of  a tapered  wing  gives  critical  values  of 
a and  d that  occur  at  divergence.  These  critical  values,  denoted  as  Ip 


I 


E 


dp  are  plotted  in  Figure  8 for  several  wing  taper  ratios.  An  approximate 
linear  relationship  again  can  be  developed  between  qQ  and  the  geometrical, 
structural  and  aerodynamic  parameters  in  the  problem.  This  relationship 
is  given  by  the  expression: 


% = 


GJr 
3 2 

aocr£  cos  A 


“ • 

JL 

er 

_Kl!l  “ Ik9  rl 


)w 


1-krtanA-Ko /GJr\  / i\  (tan  a - gf) 


(45) 


The  constants  and  are  functions  of  the  taper  ratio  a and  the  eigen- 
values found  for  each  taper  ratio.  These  constants  are  tabulated,  for 
several  values  of  a,  in  Table  1. 

TABLE  1 . - VALUES  OF  K,  AND  K?  FOR  FOUR  VALUES 
OF  TAPER  RATIO,  a = ct  /cf. 


! 

! 


A 

Ki 

K? 

0.20 

2.83 

0.614 

0.50 

2.73 

0.497 

1.00 

2.47 

0.390 

1.50 

2.22 

0.32G 

The  effect  of  taper  upon  the  divergence  of  swept  composite  wings  can 
be  seen  by  considering  the  example  for  which  all  fibers  in  the  laminate  are 
oriented  at  an  angle  e.  The  effect  on  Acr  of  the  taper  ratio  is  shown  in 
Table  2 on  page  29. 
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Figure  8 - Swept  wing_divergence  parameters  ap 
versus  dn/an  for  several  wing  taper 
ratios.  (Adapted  from  NACA  TN  No. 
1680). 
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TABLE  2.  - A VERSUS  6 FOR  THREE  TAPER  RATIOS,  X = -r- 

cr  c. 


A , critical  sweep  angle  (degrees) 

v 1 

e(deg. ) 

X=0.20 

X=0.50 

X=1 .0 

0 

1.20 

1.48 

1.89 

30 

13.71 

13.76 

13.83 

60 

35.94 

36.00 

36.08 

90 

11.13 

13.66 

17.20 

100 

-49.08 

-48.67 

-48.05 

120 

-35.45 

-35.39 

-35.30 

150 

-13.29 

-13.24 

-13.16 

The  effect  of  decreasing  taper  ratio  on  Acr  is  to  increase  Acr,  al- 
though these  increases  are  very  small.  Thus,  it  may  be  concluded  that  taper 

ratio  does  not  have  a strong  influence  on  A . 

cr 

Figure  9 illustrates  the  behavior  of  the  divergence  velocity,  (normal- 
ized with  respect  to  at  A = 0 and  e = 0),  as  a function  of  e,  for  x = 0.20. 
Comparison  of  Figure  9 with  Figure  7,  in  which  X = 1.0,  shows  that  the  sets 
of  curves  are  very  similar.  Some  differences  are  apparent  and  worth  noting. 

In  particular,  near  the  fiber  angle  9=0°,  the  divergence  velocity  at 
A = 0 declines  more  rapidly  with  e.  In  addition,  the  curve  corresponding 
to  a sweep  angle  A = -60°  does  not  attain  as  large  a maximum  value  when 
X = 0.20  as  it  does  when  x = 1.0.  It  should  be  noted,  however,  that  the 
reference  velocity  is  different  in  Figure  7 than  it  is  in  Figure  9. 
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SECTION  III  i 

CONCLUSIONS 

This  report  has  presented  a simple  theory  to  predict  the  divergence 
behavior  of  high-aspect-ratio,  laminated  composite,  swept  wings.  This 
theoretical  development  is  based  upon  laminated  beam  bending  theory  and 
aerodynamic  strip  theory.  The  equations  developed  from  this  theory  have 
been  used  to  predict  the  effect  of  laminate  construction  on  wing  divergence. 

The  conclusions  drawn  from  the  analysis  and  the  example  cases  presented 
are: 

1)  The  elastic  coupling  between  bending  and  torsion,  introduced 
by  the  laminated  material  through  the  parameter  K,  can  suc- 
cessfully negate  the  undesirable  influence  of  forward  sweep 
on  wing  divergence  for  a wide  range  of  forward  sweep  angles. 

2)  The  ratio  K/GJ  is  a very  important  parameter  to  consider 
when  designing  an  efficient  forward  swept  wing.  This  ratio 
should  be  tailored  such  that  it  has  a relatively  large  nega- 
tive value. 

3)  Wing  taper  is  important  to  the  determination  of  the  diver- 
gence speed  where  such  a speed  actually  exists.  However, 
wing  taper  is  of  secondary  importance  as  a parameter  in- 
fluencing the  design  of  a divergence  free  wing. 

4)  The  alignment  of  lamina  fibers  at  angles  from  10°  to  15° 
forward  of  the  swept  wing  box-beam  reference  axis  appears  to 
be  an  optimal  orientation  for  wing  divergence  performance. 

This  report  has  shown  that,  from  a theoretical  standpoint,  and  unlike 
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its  metallic  counterpart,  the  laminated  composite,  forward  swept  wing  is  a 
feasible  structural  design.  The  formulas  derived  and  presented  in  this 
report  provide  a quick,  inexpensive  estimation  tool  for  further  research 
on  forward  swept  wing  design. 


APPENDIX  A 


THE  DEVELOPMENT  OF  THE  STRAIN  ENERGY 
EXPRESSION  FOR  A LAMINATED  FILAMENTARY 
WING  BOX 


1.  BASIC  ASSUMPTIONS 

The  wing  structural  model  employed  throughout  this  report  is  based 
upon  the  assumption  that  the  bending  stiffness  and  torsional  stiffness 
of  the  wing  structure  are  due  entirely  to  the  presence  of  relatively 
thin,  laminated  composite  cover  sheets,  located  at  a large  distance  (in 
comparison  to  their  thickness)  from  the  wing  middle  surface.  Although 
all  of  the  bending  and  torsional  stiffness  is  assumed  to  reside  entirely 
in  the  laminated  cover  sheets,  it  would  be  a simple  task  to  add  algebra- 
ic'lly  any  additional  stiffnesses  that  might  arise  from  the  presence  of 
such  elements  as  spar  caps  or  flexible  webs. 

To  formulate  the  strain  energy  expression  that  leads  to  the  definition 
of  the  equivalent  bending  stiffness,  El,  and  torsional  stiffness,  GJ,  for 
this  type  of  box-beam  structure,  it  is  necessary  to  make  an  assumption 
about  the  deformation  behavior  of  the  wing  box.  The  well-known  Euler-Ber- 
noulli  hypothesis  that  the  strain  due  to  bending  varies  linearly  from  the 
neutral  surface  of  the  box-beam  provides  this  assumption  about  the  dis- 
placement behavior  of  the  wing.  In  addition,  it  is  further  hypothesizr-d 
that  chordwise  sections  of  the  wing,  perpendicular  to  a beam  reference 
axis,  are  rigid.  This  means  that  the  deformation  of  the  wing  box  is  a 
function  only  of  a spanwise  coordinate,  y.  In  addition  to  the  Euler-Ber- 
noulli  hypothesis,  the  material  behavior  is  assumed  to  be  linear  elastic. 


An  additional  feature  of  the  idealized  bean  ,iodel  to  be  developed  is  that 
the  reference  surface  of  the  box  beam  is  taken  at  the  geometrical  middle 
surface  of  the  wing  box.  This  convention,  commonly  used  in  laminated  plate 
theory,  is  more  convenient  than  the  alternative  of  locating  a modulus 
weighted  centroid  of  each  section  (a  modulus  weighted  centroidal  axis  will 
always  lie  in  the  neutral  surface).  If  the  laminated  box  beam  is  not  both 
elastically  and  geometrically  symmetrical,  the  middle  surface  will  not  be 
the  neutral  surface  for  bending. 

2.  COVER  SHEET  DISPLACEMENTS 

A Cartesian  coordinate  system  is  defined  in  Figure  1.  The  x-y  plane 
corresponds  to  the  geometrical  middle  surface  of  the  wing  box.  The  posi- 
tive x-ax<s  is  rearward  on  the  wing  while  the  y-axis  lies  along  the  swept 
wing  span  and  is  coincident  with  the  reference  axis  of  the  wing  box.  Dis- 
placements of  the  individual  lamina  forming  the  cover  sheets  are  given  by 
the  functions  u,  v,  w,  in  the  x,  y,  z directions,  respectively.  The  u,  v, 
w displacements  are  not  independent  functions;  they  are  functions  of  the 
variables  h and  a,  the  upward  displacement  of  the  wing  elastic  axis  (with 
respect  to  the  y-axis)  and  the  rotation  of  a chordwise  section  about  the 
y-axis,  respectively. 


In  terms  of  h(y)  and  a(y),  the  functions  u,  v,  w may  be  expressed  as: 


u(x,y)  = za(y) 

(A-l) 

v(x,y)  = vQ(y)  - z[0  - x 

(A-2?) 

or 

v(x,y)  = vQ(y)  - z[h'  - xa'] 

(A-2b) 

w(x,y)  = h(y)  - xa(y) 

(A-3) 

The  variable  vQ(y)  represents  the  spanwise  stretching  of  the  middle 
surface  caused  by  transverse  applied  loads  and  twisting. 
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3.  STRAINS  IN  THE  COVER  SHEETS 

The  strains  in  the  laminated  cover  sheets  may  be  determined  from 
Equations  (A-l)  through  (A-3).  Only  two  strains  are  nonzero;  these  strains 


eyy  = Vo  " z(h"  ■ x“"* 

Yxy  = 2za' 


(A-4) 

(A-5) 


These  strain-displacement  expressions  are  seen  to  be  consistent  with  the 
isotropic  Euler-Bernoulli  theory  of  beam  bending  in  which  there  is  neither 
transverse  shear  strain  nor  normal  stress  az  present.  The  strain-dis- 
placement expressions  are  to  be  used  to  develop  the  strain  energy  expres- 
sion for  the  box  beam,  from  which  will  result  the  El  and  GJ  expressions. 

4.  LAMINA  CONSTITUTIVE  RELATIONS 

For  a single  orthotropic  lamina,  the  relationship  between  membrane 
stress  and  strain  can  be  written,  in  matrix  notation,  as: 


^11  ^12  ^16 


°yy/  Ml2  ^22  Q26 


Q16  Q26  Q66 


(A-6) 


The  terms  Q. . are  functions  of  the  orthotrcpic  engineering  constants 

* J 

Q..  and  the  angle  e defining  the  orientation  between  the  lamina  principal 

* v' 

axes  and  the  wing  reference  axes.  The  lamina  lies  in  a plane  parallel  to 
the  x-y  plane;  the  angle  e is  as  defined  in  Figure  2.  The  terms  are 


defined  as: 


Q11  = El/^1  " u12  y21^ 
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(A-8a) 


G12  = p12  E2^  ’ p12  p21^ 


s M21  El/(1  " U21  y12^ 

(A-8b) 

Q22  * E2/(l  - y12  p21) 

(A-9) 

Q66  = G12 

(A-10) 

The  terms  Q. . can  be  defined  in  a concise  manner  [8]  by  defining  a 

* J 

; of  parameters  that  are,  in  turn,  functions  of  Q^. 

U1  = [SQ'j-j  + 3Q22  + 20^  + 4Qgg]/8 

(A-ll ) 

U2  = Wn  ■ ^22^/2 

(A-12) 

U3  = ^Qll  + G22  ' 2G12  " 4(W/8 

t A-l 3) 

U4  = ^11  + G22  + 6G12  " 4(W/8 

(A-l  4 ) 

U5  = Wll  + ^22  ’ 2Q12  + 4(W/8 

(A-15) 

terms  of  these  new  constant?,  the  matrix  terms  become 

U'n  = U-j  + U2cos2e  + U3cos4e 

(A-l 6 ) 

^12  = U4  - U3cos40 

(A-l 7) 

^22  = U1  " U2cos20  + U3cos4e 

(A-l 8) 

^16  = \ U2sin28  + U3sin40 

(A-l 9) 

^26  = 1"  u2s1n20  ' U3si’n40 

(A-20) 

^66  = U5  “ U3cos49 

(A-21 ) 

Note  that  Q-jg  and  Q2g  are  anti  symmetrical  functions  of  the  angle  0,  while 
the  other  functions  are  symmetrical. 

The  constitutive  relations  for  the  inplane  or  membrane  behavior  of 
the  individual  lamina  are  used  to  determine  the  strain  energy  functional 
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for  the  box  beam.  In  this  case,  the  individual  laminae  are  constrained  to 
act  as  a unit  because  of  the  previously  defined  strain-displacement  re- 
lations. 

5.  THE  LAMINATE  STRAIN  ENERGY  EXPRESSION 

The  strain  energy  density  per  unit  volume  for  the  box  beam  is  defined 


u*  = 2<ayy  eyy  + Txy  V 


(A-22) 


In  terms  of  this  volumetric  strain  energy  density,  a functional  U(y)  can 
be  defined.  This  latter  functional  represents  the  strain  energy  per  unit 


length  and  is  expressed  as: 


.c/2  z 


f J U*  dxdz 


(A-23) 


•c/2  -z 


The  limits  of  integration  represent:  c,  the  wing  chord;  and  zy,  the  z- 
location  of  the  upper  surface  of  the  wing  box.  Symmetry,  both  elastic  and 
geometric,  of  the  wing  box  in  the  x-direction  is  assumed  here. 

The  constitutive  relations,  Equations  (A-6),  can  be  substituted  into 
Equation  (A-23)  with  the  following  result. 


-c/2  r z 1 

U(y)  -\j  f £yy  Yxy  "^^xy'2' 


(A-24) 


-c/2  -z 


Th'  superscript  (i)  appended  to  the  lamina  elastic  constants  refers  to  those 
constants  appropriate  to  the  (i)th  layer  of  the  laminate. 

Upon  substitution  of  the  strain-displacement  expressions.  Equation 
(A-24)  becomes 


The  summation  in  Equation  (A-25)  extends  over  the  N layers  of  composite 
material.  The  constants  B-j  and  5^  are  defined  in  terms  of  lamina  coordinates 
with  respect  to  the  middle  surface.  In  terms  of  the  lamina  thickness  t^ 
and  lamina  lower  and  upper  position  coordinates  and  z^+,,  respectively, 

B.j  and  6^  are  defined  as  follows 


= J (tf  + 3tizf  + 3z1t^)  (A-26) 

and 

6,  -f  zdz 

= t.(z.+1  + Zi)/2  (A-27) 

The  term  B^  represents  the  area  moment  of  inertia  of  a strip  of  material, 
of  thickness  t^  and  unit  width,  about  the  middle  surface.  As  such,  B^  is 
always  a positive  number.  The  term  6.  represents  the  first  moment  of  the 
area  of  the  same  strip  about  the  middle  surface.  This  term  ^ is  positive 
if  the  lamina  area  centroid  lies  above  the  middle  surface  and  negative  if 
it  does  not. 

Equation  (A-25)  may  be  simplified.  It  is  to  be  noted  that  terms  con- 
taining squares  of  beam  curvature,  h",  and  the  rate  of  change  of  twist,  o', 
appear  in  this  expression.  Six  new  terms  Involving  the  laminate  elastic 
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properties  may  be  defined,  leading  to  a more  concise  definition  of  Equa- 
tion (A-25).  These  terms  are: 


(A-28) 


(A-29) 


(A-30) 


(A-31 ) 


(A-32) 


(A-33) 


(A-34) 


With  these  newly  defined  constants,  U(y)  may  be  written  as: 


_r  SL  fc-55  rs®**?  *Zrf  ^ 


w *?^v  -«??^rrsa^5 


Equation  (A-35)  reveals  that  bending  deformation  and  torsional  de- 
formation are  coupled  through  the  constant  KQ»  defined  in  Equation  (A-31), 
and  that  additional  coupling  between  the  middle  surface  stretching  and  both 

bending  and  torsion  also  occurs.  The  terms  El  and  GJ  are  defined  as 

oo 

they  are  because  they  occur  in  the  energy  terms  associated  with  uncoupled 
bending  and  torsion,  respectively. 

In  Reference  7.  the  term  in  the  strain  energy  expression  containing 
SQ  is  shown  to  be  insignificant  in  comparison  to  that  term  involving  the 
term  GJQ  as  long  as  the  wing  has  a moderate  to  high  aspect  ratio.  This 
term  is  also  ignored  in  the  analysis  in  the  present  study. 
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APPENDIX  B 

THE  DEVELOPMENT  OF  THE  COMPOSITE  BEAM  BENDING 
AND  TWISTING  EQUILIBRIUM  EQUATIONS  FOR  DIVERGENCE 

1.  GOVERNING  EQUATIONS  FOR  A BEAM  WITH  CONSTANT  CROSS-SECTION 

The  strain  energy  functional  can  be  used,  in  conjunction  with  the 
Principle  of  Virtual  Work,  to  derive  the  equilibrium  equations  necessary 
to  solve  the  divergence  problem.  For  the  present  problem,  the  Principle 
of  Virtual  Work  can  be  stated  as  follows: 

<$  / U(y)dy  = / (py5vQ  + t 5a  + p5h)dy  (B-l) 

*4  *4 

The  terms  py,  t and  p in  Equation  (B-l)  represent  the  applied  generalized 
forces  per  unit  length  associated  with  the  degrees  of  freedom  vQ,  a,  h, 
respectively.  These  generalized  forces  are  distributed  along  the  length 
of  the  beam;  no  concentrated  forces  or  moments  are  applied.  For  our  analy- 
sis, the  generalized  force  Py  is  zero,  since  it  represents  an  external 
axial  load  term. 

Manipulation  of  Equation  (B-l)  yields  three  differential  equations  of 
equilibrium  (in  terms  of  displacement  derivatives)  together  with  the  neces- 
sary boundary  conditions.  These  differential  equations  are,  for  the  con- 


stant  property  beam  model , 

as  follows: 

'V" 

-A  v " + R h'" 
22  0 22 

= 0 

(B-2) 

• Vo"  + V" 

- t(y) 

(B-3) 

- V"  - Vo"' 

= p(y) 

(B-4) 

The  boundary  conditions  at 

y * 0 are  as  follows: 

vQ(0)  = 0 
h(0)  = h* (0)  = 0 

a(0)  = 0 

At  y = a,  the  boundary  conditions  are: 


6V  + B33Vo'  - Koh"  = 0 

B33a'  + A22Vo'  " B22h"  = 0 
EI0h"  - K0a‘  - B22V0'  = 0 

EIoh’"  - V"  - B22vo"  = 0 


(B-5 ) 

(B-6,7) 

(B-8) 

(B-9) 
(B-10) 
(B-ll ) 
(B-12) 


An  understanding  of  the  physical  significance  of  the  equations  of 
equilibrium  can  be  found  by  formulating  expressions  for  internal  force  and 
moment  resultants  in  terms  of  the  displacement  derivatives.  Let  Py  repre- 
sent the  internal  resultant  axial  force  present  in  the  beam.  This  force 
resultant  is  defined  as: 


If 


py  JJ  °yydxdz  = A22Vo'  ' B22h"  + B33“' 

Area 

The  internal  bending  moment  Mx  is  defined  as: 


(B— 13) 


V 


-jr 


vdxdz 


* -B-,v  ' + El  h“  - K o' 
22  0 0 0 


(B-14) 

From  static  equilibrium  of  an  infinitesimal  beam  section,  the  internal  shear 


force  resultant  Vz  is  given  by  the  equation 

dM 

Vz  " ~i  - "B22vo"  + EIoh'"  * Ko“" 
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(B-15) 
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Finally,  the  internal  torque  T is  defined  as: 


T = 


// 

Area 


(2Txy  ' *TxZ>dxdz 


(B-16) 


= GJ  o'  + B,,v  ' - K h" 
0 33  0 0 


Comparison  of  Equation  (B-13)  with  Equation  (B-2)  shows  that  the 
latter  equation  may  be  written  as: 


dPy 

-/■  = 0 
dy 


(B-17) 


Since  there  is  no  loading  in  the  axial  direction,  Equation  (B-17)  can 
be  integrated  once  to  produce  the  following  result. 


Py  = 0 = A22vo ‘ ' B22h"  + B33a’ 


(B-18) 


Equation  (B-18)  may  be  solved  for  vQ' 


B22  B33 

\/  1 - LCm  h11  ^ * 

° '*22  " hi 


(B-19) 


Equation  (B-19)  may  be  differentiated  to  obtain  identities  for  higher  order 
derivatives  of  v , h and  a.  The  identities  obtained  by  differentiating 

J 

Equation  (B-19)  may  be  substituted  into  the  equation  for  bending.  Equation 
(B-4)  and  the  equation  for  torsion.  Equation  (B-3).  The  following  equations 
result. 


EIo- 


<h£ 

A22 


Ko  " 


B22  B33 
A22 


jv 


h'"  + 


Ko  ' 


B22  B33 
A22 


= p(y) 


(s33y 

GO  - - v-  ■ 
o A22 


a"  = 


-t(y) 


(B-20) 


(B-21 ) 


These  equations  are  independent  of  the  middle  surface  stretching  term  vQ. 
The  boundary  conditions  for  the  problem  may  also  be  modified  so  that 
derivatives  of  vQ  are  eliminated. 

Definitions  of  reduced  bending  stiffness,  El,  torsional  stiffness,  GI, 


and  the  coupling  parameter  K are  readily  available  from  Equations  (B-20) 
and  (B-21 ) . With  these  definitions.  Equations  (B-20)  and  (B-21)  become: 


El  h1V  - K a"'  = p(y) 

(B-22) 

-K  h"'  + GJa"  = -t(y) 

Boundary  conditions  for  these  equations  are  as  follows: 

(B-23) 

at  y = 0,  a(0)  = 0 

(B-24) 

h ( 0)  = h ' (0)  = 0 

(B-25) 

at  y = £,  GJa1  - K h"  = T(£)  = 0 

(B-26) 

El  h"  - K a'  = Mx(£)  = 0 

(B-27) 

El  h'"  - K a"  = Vz(£)  = 0 

(B-28) 

2.  GOVERNING  EQUATIONS  FOR  A BEAM  WITH  A NON-UNIFORM  CROSS-SECTION 

The  Principle  of  Virtual  Work  also  can  be  used  to  derive  the  governing 
equations  for  a nonuniform  beam.  However,  an  alternative  way  to  derive 
these  equations  is  to  use  the  definitions  of  the  bending  moment,  twisting 
moment  and  resultant  shear  force,  in  terms  of  the  deformations  h,  a and 
vQ,  derived  in  the  previous  section. 

Onc^  .again  the  axial  force  equation  yields  a relationship  identical  to 
that  given  in  Equation  (B-19).  The  bending  equation  reads 

H22V  + EIoh"  ' Koa,)"  = p(y)  (B_29) 

while  the  torsion  equation  becomes 

( B33vo ' ’ Koh"  + ™oa'y  = 't(y)  (B'30) 
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Elimination  of  the  axial  stretching  term,  vQ,  results  in  two  equations 
involving  only  two  variables,  h and  a.  These  equations  are: 

[El  h"  - K a1]"  = p(y)  (B-31) 

[-K  h"  + GJ  «•]'  = -t(y)  (B-32) 

The  definitions  of  the  terms  El,  GJ,  and  K are  identical  to  those  discussed 
in  the  previous  section.  Boundary  conditions  for  Equations  (B-31)  and 
(B-32)  are: 


at  y = 0 o(0)  = h(0)  = h’{0)  = 0 

at  y = l Elh"  - K o'  * 0 * M (A) 

[El  h"  - K o’]*  = 0 = V (A) 
GJ  a'  - K h,!  = 0 = T(a) 

3.  AIRLOADS  FOR  THE  DIVERGENCE  PROBLEM 


(B-33,34,35) 

(B-36) 

(B-37) 

(B-38) 


Figure  3 illustrates  the  geometry  employed  to  determine  the  incremental 
lift  forces  and  pitching  moments  on  the  swept  wing.  The  airloads  are  modelled 
from  strip  theory  assumptions,  assumptions  that  tend  to  overestimate  the 
aerodynamic  loading  and  also  give  values  for  the  airloads  that  are  somewhat 
inaccurate  at  large  sweep  angles.  Nevertheless,  the  value  of  strip  theory 
aerodynamics  lies  in  its  ability  to  provide  reasonably  accurate  aeroelastic 
trend  information,  allowing  the  establishment  of  the  proper  relationships 
between  the  geometric,  aerodynamic  and  structural  parameters  in  an  aero- 
elastic problem. 

References  1 and  9 detail  the  development  of  the  use  of  strip. theory 
aerodynamics  for  swept  wing  aeroelasticity  problems.  Of  importance  to  the 
divergence  problem  is  the  determination  of  the  incremental  lift  per  .unit 
of  span  (along  the  swept  y-axis)  arising  from  the  deformations h and  a.  The 
reader  is  reminded  that,  for  this  study,  h and  o represent  the  perturbation 
deformations  away  from  some  slightly  deformed  static  equilibrium  position. 


\ 


\ 

\ 
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The  lift  per  unit  span  due  to  these  deformations  is  given  by  the  expression 
p{y)  = {caQqcos2A)  (a  - h'tanA)  (B-39) 

While  the  pitching  moment  about  the  y-axis  is  written  as: 

t(y)  = (qcaQe  cos2a)  (a  - h’tanA)  (B-40) 

where  aQ  = lift  curve  slope  (two-dimensional). 

c = chordwise  dimension  perpendicular  to  the  y-axis. 
e = the  distance  between  the  line  of  aerodynamic  centers  (A.C.) 

and  the  y-axisv  measured  positive  aft  from  the  A.C. 

A = the  sweep  angle  of  the  reference  axis,  positive  rearward. 

4.  NONDIMENSIONALIZED  GOVERNING  EQUATIONS  FOR  A UNIFORM  WING 

With  the  inclusion  of  the  airloads  given  in  Section  B3,  the  coupled 
equations  governing  divergence  of  a uniform  cross-section,  composite  swept 
wing  are: 

El  hiv  - K a'"  = (caQqcos2A)(a  - h'tanA)  (B-41) 

and 

GJ  a"  - K h'"  = -(qca0ecos2A)(a  - h'tanA)  (B-42) 

These  equations  may  be  written  conveniently  in  nondimensional  form.  To 
accomplish  this,  new  variables  are  defined  as  follows: 

n = y/£  (B-43) 

w = h/£  (B-44) 

r = dw/dn  = dh/dy  (B-45) 

These  definitions  lead  to  the  following  nondimensional,  coupled  equilibrium 
equations. 

r'"  - k a"'  * a(a  - rtanA)  (B-46) 

a"  - gr"  = -b(a  - rtanA)  (B-47) 

46 


where  ( )'  = d{  )/d0.  The  boundary  conditions  for  these  equations  are 

written  as  follows: 

r(0)  = 0 

(B-48) 

P 

O 

It 

O 

(B-49) 

r'(l)  - kct'(l)  = 0 

(B-50) 

a * (l ) - gr'(l)  = 0 

(B-51 ) 

r"(l)  - ka"(l)  = 0 

(B-52) 

The  nondimensional  parameters  that  appear  in 

Equations  (B-46) through  (B-52) 

are  defined  as  follows: 

3 2 

a = qc£  aQcos  a/EI 

(B-53) 

2 2 

b = qcji  aQecos  a/GJ 

(B-54) 

k = K/EI 

(B-55) 

g = K/GJ 

(B-56) 

5.  NONDIMENSIONAL  EQUILIBRIUM  EQUATIONS  FOR  A TAPERED  WING  PLANFORM 

A wing  whose  chord  is  linearly  tapered  along  the  wing  span  and  whose 
cross-sectional  dimensions  are  proportional  to  the  local  chord  will  have 
bending  stiffness  and  torsional  stiffness  distributions  that  vary  as  the 
fourth  power  of  the  chord;  the  coupling  parameter  K will  also  vary  as  the 
fourth  power  of  the  chord.  To  begin  the  development  of  the  equilibrium 
equations  applicable  to  this  problem,  let  us  define  the  following  terms: 


c = fc * = fc 

root  r 

(B-57) 

where 

f = 1 - n(l  - A) 

(B-58) 

and 

A = ctip/croot  = ct/cr  = wing  taper  ratio 

(B-59) 

The  subscript  notation  ( )r  refers  to  values  evaluated  at  the  wing  root, 
while  n = y/i. 
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Since  the  chord  varies  linearly  with  n and  all  wing  sections  are  geo- 
metrically similar,  the  following  stiffness  parameter  relationships  result. 


El  = f*EI 


GJ  = rGJ. 


(B-60) 

(B-61) 

(B-62) 


Equations  (B-31 ) and  (B-32)  govern  the  static  equilibrium  of  the  nonuniform 
wing  and  can  be  combined  with  Equations  (B-39)  and  (B-40)  (note  that,  in 
the  latter  two  equations,  c = fcj  1 to  yield  an  equation  that  is  a function  of 
the  nondimensional  independent  variable  n.  The  variable  f = 1 - n(l  - a) 
also  can  be  chosen  as  the  independent  variable.  Since  the  latter  choice  of 
an  independent  offers  substantia!  mathematical  advantages  when  attempting 
a closed  form  solution  to  the  problem,  f,  not  is  chosen  as  the  independent 
variable.  With  f as  the  independent  variable,  the  bending  equation  reads 
(note  that  ( )'  = d(  )/df  below): 

f3[r"'  - kra'"]  + 8f2[r"  - kra"]  + 12f[r'  - kra']  ‘ V + taHX  s 0 (B"63) 


where 


qcr*  aQsinAcosA 


EIrO  - A}; 


(B-64) 


The  torsion  equation  becomes 


f2[a”  - grr"]  + 4f[a'  - grr']  + aT[u  - rtanA]  - G 

where  2 ? 

qcrerra0cos^A 

aT  ‘ ur(l  - X)2  " 

Boundary  conditions  at  f = 1 (n  = 0)  are, 

a(l)  = o 
r(l)  » 0 
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(B-65) 

(B-66) 

(B-67) 

(B-68) 


while,  at  f = A (n  = 1), 

M(A)  = r ‘ (X)  - kfa'(A)  = 0 (B-69) 

V( A)  = 4A3( r 1 (A)  - kra'(A))  + A4(r"(A)  - kron(A))  = 0 (B-70) 

T(a)  = o' (a)  - gr * ( A ) = 0 (B-71) 

While  these  latter  equations  and  boundary  conditions  now  appear  rather 
complex,  they  can  be  reduced  to  a form  in  which  they  can  be  readily  solved. 
This  simplification  and  solution  is  discussed  in  Appendix  C. 
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SOLUTION  OF  THE  WING  DIVERGENCE  EQUATIONS  FOR  CONSTANT 


CHORD  AND  LINEARLY  TAPERED  WINGS 


1 . SOLUTION  FOR  CONSTANT  CHORD  WINGS 


In  Section  B4,  two  coupled  differential  equations,  Equations  (B-46) 


and  (B-47),  were  developed;  these  equations  govern  the  divergence  of  a 


swept  wing  with  a constant  chord.  The  solution  of  these  equations,  sub- 


ject to  boundary  conditions.  Equations  (B-48)  through  (B-52),  is  detailed 


in  this  appendix. 


Let  us  begin  by  defining  two  new  dependent  variables,  <j>(n)  and  V»(n) , 


to  be  used  to  transform  Equations  (B-46)  and  (B-47)  into  a single  equation 


in  terms  of  the  variable  ag  = a - r tanA.  In  terms  of  $ and  *,  the  variables 


r and  a are  defined  as: 


r = 4>  (n ) + k*(n) 


(C-l) 


a = if/(n)  + g*(n) 


(C-2) 


where  the  constants  k and  g have  been  previously  defined.  Substitution  of 


Equations  (C-l)  and  (C-2)  into  Equations  (B-46)  and  (B-47)  yields  the  fol- 


lowing: 


♦ - [r.-akg][(l  - k tanA)*  + (g  - tanA)*]  = 0 


(C-3) 


* " + []  kg][(l-  k tanA)*  + (g  - tanA)*]  = 0 


(C-4) 


Substitution  of  expressions  for  *(r,)  and  *(n),  as  functions  of  r and  o. 


obtained  from  inverting  Equations  (C-l)  and  (C-2),  into  the  expression 


contained  within  the  second  set  of  brackets  in  both  Equations  (C-3)  and 
(C-4)  gives  the  following  result. 


Mm 


(1  - k tanA)i|i  + (g  - tanA)«j>  = a -rtanA  (C-5) 

The  variable,  ag,  is  defined  as 

ag  = a - r tanA  (C-6) 

The  term  afi  is  the  local  wing  angle  of  attack,  perpendicular  to  the  elastic 
axis,  due  to  wing  bending  and  torsional  deformation.  Equations  (C-3)  and 
(C-4)  now  read 

- [Tr^E8T]  = 0 <«> 

♦'  + [(1  -b|<g)]  [ae]  ‘ 0 (C'8) 

If  we  first  differentiate  Equation  (C-8),  multiply  it  by  the  factor 
(1  - k tanA)  and  then  add  the  result  to  the  product  of  Equation  (C-7) 
times  the  factor  (g  - tanA),  the  result  is  a single  equation,  written  as: 


Equation  (C-9)  is  a third  order  differential  equation  with  constant  coef- 
ficients. 

Boundary  conditions,  in  terms  of  the  variable  ag,  are  constructed  in 
the  following  manner.  Equations  (B-48)  and  (B-49)  can  be  combined  to  give 
the  equation 

a(0)  - r(0)  tanA  = ag(0)  = 0 (C-10) 

An  examination  of  Equations  (B-50)  and  (B-51)  shows  that  both  a'(l)  and 
r!(l)  must  oe  zero,  thus,  the  following  boundary  condition  results. 

o'(l)  - r'(l)  tanA  = <^{1)  = 0 (C-11) 

Next,  Equation  (B-47)  is  evaluated  at  n = 1»  multiplied  by  the  factor 
(1  - k tanA),  and  then  added  to  the  result  of  the  product  of  Equation  (B-52) 
and  the  factor  (g  - tanA).  This  operation  provides  the  third  boundary  con- 
dition, Equation  (C-12). 
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Equations  (C-9)  through  (C-12)  constitute  an  eigenvalue  problem.  The  solu- 


tion to  Equation  (C-9)  is  of  the  form 


a.  = Ce 
e 


(C-13) 


Substitution  of  Equation  (C-13)  into  Equation  (C-9)  yields  a characteristic 
equation  whose  solution  provides  three  values  for  the  coefficient  r.  This 


equation  is 


r^  + ar  + d - 0 


where  the  constants  a and  d are  defined  as: 

r-  h - k tantl  fqcet2  aocos2j 
a ' n GJ 


- _ tanA  - < 
d ' T^kg 


llf^  3 

eT 


3 2 

aQcos  A 


(C-14) 


(C-15) 


(C-16) 


Once  the  three  roots  of  Equation  (C-14)  are  determined  as  functions 
of  a and  d,  Equation  (C-13)  may  be  used  to  construct  the  stability  determi- 
nant stemming  from  enforcement  of  the  homogenous  boundary  conditions  for 
the  problem.  The  solution  of  the  stability  problem  is  accomplished  by 
finding  values  of  a and  d which  render  this  stability  determinant  zero.  These 
critical  values  of  a"  and  d are  denoted  as  a^  and  d^,  respectively. 

The  numerical  solution  to  the  stability  problem  itself  is  lengthy.  For- 
tunately, Reference  1 details  the  solution  to  an  identic’  set  of  equations. 
Although  the  problem  in  Reference  1 involves  a metallic  wing,  the  solution 
presented  in  Reference  1 can  be  adapted  to  the  composite  wing  by  a suitable 
change  in  the  definition  of  appropriate  constants.  An  approximate  solution 


is  given  in  Reference  1 and  discussed  in  Section  2.4  of  the  present  report. 


This  approximate  solution  Is  as  follows: 

aD  = 2.47  + 0.390dD  (C-17) 

From  Equation  (C-17),  the  value  of  q^,  the  dynamic  pressure  at  wing 
divergence,  is  found  to  be 


% = 


fiii?  (i  - Ml 

El 

3^  2 

c£  aQcos  A 

1. -ktanA  n 

/GO,  0.39 (tanA  - g) 

LV  *EI;  J 

(C-18) 


2.  SOLUTION  FOR  LINEARLY  TAPERED  WINGS 

In  Section  5 of  Appendix  B,  the  governing  equations  for  a wing  with 
linear  taper  are  presented.  The  solution  to  these  equations  is  detailed 
in  this  section.  Equations  (B-63)  and  (B-65)  may  be  combined  in  a manner 
similar  to  that  outlined  in  Section  1 of  Appendix  C.  The  procedure  is  as 
follows.  First,  the  two  variables  $ and ^ are  defined.  These  new  variables 
are  defined  as  follows: 

r = + krij>(n)  (C-19) 

a = 4>(  n)  + 9r<j>(n)  (C-20) 

In  addition,  we  define,  as  in  the  case  of  the  constant  chord  wing, 

ae  = a - r tanA  (C-21) 

with  the  definitions  contained  in  Equations  (C-19)  through  (C-21)  and  with 
( )'  = d(  )/df,  Equation  (B-63)  becomes 

fV"  * 8fV  ♦ 12V  * (1  ■ 0 (C-22) 

Similarly,  Equation  (B-65)  becomes 

f V + 4f*'  + g = 0 (C-23) 

Kr9r 
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The  following  procedure  results  in  the  reduction  of  Equations  (C-22)  and 
(C-23)  to  a single  equation  in  terms  of  og: 

1)  Multiply  Equation  (C-22)  by  (gr  - tanA). 

2)  Differentiate  Equation  (C-23)  once  with  respect  to  f and 
multiply  the  result  by  the  quantity  f(l  - krtanA). 

3)  Add  the  results  from  (1)  and  (2)  together. 

4)  Make  use  of  identities  to  reduce  the  equation  resulting  from 
(3)  to  the  form 

f 3(31  g + Sf2^  + (12  + aT)fa;  + (2?t  - dT)ae  = 0 (C-24) 

where 


1 - k tanA 
r 

O A 

qercra  aocos  A] 

1 - kg 
rJr 

- GJr(l  - A)2  -I 

tanA  - g„ 
r 

3 ? 

qc„JT  a cos  A 
^ r o 

1 ' Mr 
r r 

EIr(l  - X)3  . 

(C-25) 


(C-26) 


The  boundary  condition  at  f = 1 (the  wing  root)  is: 

ae(l)=0  (C-27) 

while,  at  f = X,  the  wing  tip,  the  bending  moment  and  torque  conditions, 

M(x)  = 0 and  T(x)  = 0,  give 

a;(x)  = 0 (C-28) 

The  condition  of  zero  shear  at  the  wing  tip,  when  combined  with  Equa- 
tion (B-65),  evaluated  at  f * X,  yields 

X^ct”  ( X ) + ?Toe(x)  = 0 (C-29) 

The  solution  to  Equation  (C-24)  subject  to  the  boundary  conditions  out- 
lined above  follows  that  described  in  Reference  1 for  a metallic  tapered  wing. 
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Reference  1 presents  an  approximate  linear  relationship  between  the 
critical  value  of  dynamic  pressure,  qD>  and  the  geometric,  aerodynamic  and 
structural  parameters  in  the  tapered  wing  problem.  This  solution,  when 
modified  to  fit  the  composite  wing  problem,  provides  the  following  equation. 


“r  1 

/ l \ 

^(1  - krgr) 

qD  = 

3 — T 

a c £ cos  A 
or 

\er  ) 

.1  -krtanA-K2(^)(i)(tanA-gr). 

The  constants  K^  and  Kg  are  functions  of  the  wing  taper  ratio  *.  These 
constants  were  calculated  for  this  study  and  are  tabulated  below  for  four 
taper  ratios.  The  values  given  here  differ  slightly  from  those  found  in 
Reference  1.  The  reason  for  this  slight  difference  is  probably  due  to  the 
improved  accuracy  available  today  on  modern  computers. 


\ 

Ki 

0.20 

2.83 

0.614 

0.50 

2.73 

0.497 

1.00 

2.47 

0.390 

1.50 

2.22 

0.326 

3.  CORRECTION  OF  AERODYNAMIC  LOADS  FOR  ASPECT  RATIO  AND  SWEEP 

In  Reference  1,  a correction  for  the  effects  of  aspect  ratio  and  sweep 
is  suggested.  This  formula  Is 

Sa  = a0W(flR  + 4 cos  (C-31) 

where  aQ  is  the  2-D  lift  coefficient  for  the  unswept  wing  section.  This 
correction  may  be  used  in  Equation  (C-18)  and  in  Equation  (C-30)  in  place 
of  aQ.  This  correction  is  used  throughout  this  report  for  all  examples 
presented. 
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